ARCHIMEDEAN-TYPE FORCE IN A COSMIC DARK FLUID: 
I. EXACT SOLUTIONS FOR THE LATE-TIME ACCELERATED EXPANSION 
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We establish a new self-consistent model in order to explain from a unified viewpoint two key fea- 
tures of the cosmological evolution: the inflation in the early Universe and the late-time accelerated 
expansion. The key element of this new model is the Archimedean-type coupling of the dark matter 
with dark energy, which form the so-called cosmic dark fluid. We suppose that dark matter particles 
immersed into the dark energy reservoir are affected by the force proportional to the four-gradient 
of the dark energy pressure. The Archimedean- type coupling is shown to play a role of effective 
energy-momentum redistributor between the dark matter and the dark energy components of the 
dark fluid, thus providing the Universe evolution to be a quasiperiodic and/or multistage process. 
In the first part of the work we discuss a theoretical base and new exact solutions of the model 
master equations. Special attention is focused on the exact solutions, for which the scale factor is 
presented by the anti-Gaussian function: these solutions describe the late-time acceleration and are 
characterized by a nonsingular behavior in the early Universe. The second part contains qualita- 
tive and numerical analysis of the master equations; we focus there on the solutions describing a 
multi-inflationary Universe. 
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I. INTRODUCTION 

The concepts of dark energy (DE) and dark matter 
(DM) (see, e.g., [l|-[3{ and [JQ for review and refer- 
ences) are the basic elements of modern cosmology and 
astrophysics. These elements were introduced into the 
scientific lexicon in two different ways: dark energy is 
considered to be a reason for the late-time accelerated 
expansion of the Universe 0, H[ , while dark matter is usu- 
ally associated with the explanation of the fiat velocity 
curves of the spiral galaxies rotation [l,[l(|. Nevertheless, 
there exists a tendency to consider dark energy and dark 
matter as two manifestations of one unified dark fluid 
(see, e.g., [TT| - [l6| ). The models of interaction between 
two constituents of the dark fluid, as well as the models 
of interactions of dark energy and/or dark matter with 
the standard (baryon) matter, are the subjects of wide 
discussion. The total contribution of dark energy and 
dark matter into the energy balance of the Universe is 
estimated to be about 95%. Thus, the coupling between 
these two constituents of the dark fluid seems to be the 
most important element in the list of cosmic medium 
interactions, and the dark fluid can be considered as a 
thermodynamic reservoir for baryon matter. 

The most developed model of the coupling between 
DM and DE constituents of the dark fluid is based on the 
two- fluid representation of the cosmic medium (see, e.g., 
[l7l - |22j ] ) . In this approach the interaction terms ±Q ap- 
pear in the right-hand sides of separate balance equations 
for the DE and DM with opposite signs and disappear in 
a sum, when one deals with the total balance equation. 
The modeling of the coupling term Q has in most cases 
a phenomenological character and is based on the ansatz 
that Q is a function (e.g., linear or power- law) of the en- 



ergy densities of the DE and DM, of the Hubble function, 
H, etc. We formulate the theory of interaction between 
DE and DM using the relativistic hydrodynamics for dark 
energy and relativistic kinetics for dark matter. We sug- 
gest that the DE acts on the DM particles by means of 
some effective force, and the corresponding model force 
four-vector is introduced into the kinetic equation. The 
backreaction of the DM on the DE is described using 
a force-type term in the hydrodynamic equations. The 
total system (DM plus DE) is considered to be conserva- 
tive. The concept of the Archimedean-type force can be 
naturally generalized for the description of the DE action 
on the baryon matter; however, now we restrict ourselves 
by the model of the DE and DM interaction. 

We discussed the structure and properties of various 
effective forces, which appeared in the cosmological con- 
texts, in the papers [23^-28], the relativistic generaliza- 
tions of the Stokes force, Langevin force, antifriction and 
tidal forces being investigated in detail. Concerning the 
force acting on the DM particles from the DE we would 
like to introduce the so-called Archimedean-type force. 
This choice can be motivated as follows. First of all, 
this force is a relativistic generalization of the classical 
Archimedean force, proportional to the three-gradient of 
the Pascal pressure; thus, the model under discussion 
is based on the well-known and well examined scheme 
of interaction. Second, the DE pressure is assumed to 
be of the same order as the DE energy density (73% 
of the Universe energy density), thus, the Archimedean 
effect on the dark matter could be significant. Third, as- 
suming that the DE pressure can be negative, we obtain 
that such a force can be attractive in contrast to ex- 
pulsive classical Archimedean force. Admitting that the 
DE pressure changes the sign in the course of the Uni- 
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verse evolution, one can describe a multistage (or even 
quasiperiodic) character of the cosmological expansion, 
for which epochs of deceleration are changed by epochs 
of acceleration and vice versa. (The interest in models of 
this type was renewed by the paper (29l|V 

We show that, in principle, the Archimedean-type force 
can effectively redistribute 95% of the Universe's energy 
between the DE and DM constituents, thus guiding the 
time evolution of the cosmic medium as a whole. We 
divided the work into two parts: the first one contains 
pure analytical results and exact solutions of the model; 
in the second part we focus on the numerical and quali- 
tative analysis of the model. 

The first part of the work is organized as follows. 
In Sec. II we derive the master equations of the model 
with an Archimedean-type force. In particular, in Sec. 
II. B based on the kinetic approach we introduce the 
Archimedean-type force, obtain basic integrals of mo- 
tion, construct the distribution functions and calculate 
their macroscopic moments as functions of the DE pres- 
sure. In Sec. II. C we formulate the balance equation 
for the dark fluid. In Sec. II. D we discuss the extended 
(inhomogeneous) equation of state for the dark energy 
and introduce the key equation for the DE pressure evo- 
lution. Sec. Ill contains discussions about two classes 
of exact solutions. In Sec. III. A we consider a special 
(constant) exact solution for the case when the guiding 
parameter of the model a is not equal to its critical value, 
i.e. a 7^ — 1, and we analyze the problem of asymptotic 
stability of this solution. In Sec. III.B we focus on the 
special case a = — 1 and obtain exact solutions of the 
anti-Gaussian type for the following submodels: (i) the 
massless DM, (ii) the cold dark matter, and (hi) the sub- 
model with the DE domination. In Sec. IV we discuss 
obtained analytical results. 



II. MASTER EQUATIONS 

We consider a cosmic medium, which consists of two 
interacting components. One of them, the dark matter, 
can be described in the framework of general relativistic 
kinetic theory [3(| Hl|. The second component, dark en- 
ergy, is considered as a perfect fluid with inhomogeneous 
equation of state. These two components interact grav- 
itationally, i.e., both of them contribute the energy and 
momentum terms into the total stress-energy tensor of 
the system as a whole, indicated as the dark fluid. In 
addition, we assume that DE and DM interact by means 
of force of the Archimedean type. Mathematically, this 
model can be described as follows. 



A. Equations for gravity field 



model with the metric 

ds 2 = dt 2 - a 2 {t)[{dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 } . (1) 
Here and below c = 1. The Einstein equations 
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Rik - ^Rgik = K T^ total) + Ag lk , (2) 



can be reduced to the well-known system 



a 47rG r „ A 

- = — 5-P + 3P+-T. 
a 6 3 



irG „ A , . 

£ + ( 3 ) 



a J 3 3 



where the dot denotes the derivative with respect to time, 
£{t) is the total energy of the system as a whole, and 
V is the corresponding pressure. Since the cosmological 
constant A is frequently interpreted in terms of vacuum 
energy and can be considered as a dark energy candidate, 
we assume the following decompositions of the total en- 
ergy and pressure 



A 
8^G 



= P + E, V- 



A 
8^G 



= U + P. 



(4) 



In such decomposition the functions p(t) and Tl(t) de- 
scribe the dark energy, while the functions E(t) and P(t) 
describe the dark matter; thus, the term -J^i is consid- 
ered to be incorporated into the DE energy density p and 
DE pressure II. Now the master equations for gravity 
field are 



a 
a 



4ttG 



[(p + E)+3(U + P)} 



tt9 8ttG , 

H 2 = —(p + E), 



(5) 



(6) 



where, as usual, H(t) = - is the Hubble function. 



B. Kinetic equation for relativistic DM particles 

Let us consider the general relativistic kinetic equa- 
tions 



'(a) 



d 

dx 1 



d 



for the distribution functions // a ), describing the evolu- 
tion of the DM particles with the masses m( a ) . It is worth 
mentioning that the DM can consist of a few sorts of par- 
ticles (massless and massive), thus, we use the index (a) 
to distinguish them. The kinetic equations ([7]) are of the 
collisionless type, i.e., we neglect direct interparticle col- 
lisions, but consider the force four-vector J 7 ^ to guide 
the DM particle dynamics. Two subsets of characteristic 
equations 



We consider the spatially homogeneous Friedmann- 
Lemaitre-Robertson- Walker (FLRW) cosmological 



dp 



:dx 



ds 



ds 



dx 1 
ds 



Ha) 



(8) 
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show that the DM particle motion is not geodesic because where p° is the DM particle energy 
of the force which appears as a result of the DE 

action. The third subset of characteristic equations 



Po = \ m 



df( a 



ds 



dp k 



(9) 



demonstrates the important role of the divergence dpk 
in the evolution of the distribution functions. 



(a) 



V Vo 



(13) 



(Greek indices run from 1 to 3). Thus, the Archimedean- 
type force (fTTj) is divergence-free, if and only if the DE 
pressure IT is constant. 



2. Integrals of motion 



1. Archimedean-type force 

In classical physics the term "Archimedean force" ap- 
pears when one deals with a body immersed into a nonho- 
mogeneous liquid. This three-dimensional force -F(Arch) 
is proportional to the spatial gradient of the Pascal pres- 
sure P( P 

ascal) 



^(Arch) — "^)VP(p asca l) 



(10) 



Clearly, in the Archimedean historical experiment the 
gradient of the pressure was equal to the product of the 
water mass density P( wa t cr ) and the free-fall acceleration 
g in Syracuse. The coefficient Vq relates to a body vol- 
ume and can be represented as a quotient 



Ma 



, where 

P(body) 

-^(body) is a mass of the body and P(body) is its mass 
density. The minus sign corresponds to the fact that the 
Archimedean force is the buoyancy (expulsive) one. 

This classical force can be generalized for the case of 
rclativistic DM particle immersed into the DE reservoir. 
The force four- vector 



Using the force four-vector given by (jlip one can ex- 
tract the following self-closed subsystem from the char- 
acteristic equations (JSJ) : 



dt 



-V (a) p Q n . 



dt 
ds 



'(a) 



(14) 



These equations yield immediately three first integrals 
(for the covariant components of the particle momentum) 

Pa (t)=C a exp {V (a) p(i )-n(t)]} , Pa (t )=C a . (15) 

Thus, any deviation of the DE pressure Tl(t) from 
the initial value II(to) generates particle accelera- 
tion/deceleration due to the Archimedean- type interac- 
tion. The energy of the DM particle of the sort (a) 




a(tp) 
a(t) 



o 2V (a) [n(t )-n(i)] 



(16) 



deviates from its initial value 



^(a) = TO(a)V (a) 



p"p n 



p l pi 



(ii) 



can be considered as a generalization. First of all, it is 
a force proportional again to the gradient of pressure II, 
but now we deal with a four-gradient instead of gradient 
three-vector in classical physics, and IT is the DE pres- 
sure instead of the Pascal one. Second, the projector in 
the square brackets guarantees that the force four- vector 
is orthogonal to the particle momentum four- vector, i.e., 
F l ^Pi=Q. The last property provides the existence of the 

first integral of motion p k pk=m 2 a ^=const, which guaran- 
tees the particle mass conservation. New constants V( a ) 
have the dimensionality of inverse energy density (re- 
member that c=l) and describe some new effective con- 
stants of interaction. The sign in front of these constants 
is positive, and this is in agreement with the definition 

(ITU)) and the accepted signature (H ) of the metric 

(jTJ) . This force does not belong to the class of gyroscopic 

<9.F fc 

forces, for which g p V =0- Indeed, simple calculations 



give 



8F k 
dp k 



= -3m (a) V( a ) 



P fc v fc n 
p 1 pi 



= -3V, 



(a) 



P °ii 



in 



(a) 



7^0, (12) 



P°(to) 



where 



m 2 a) a 2 (t ] 



[Cl 



cl 



CI 



(17) 



(18) 



due to the following effects. The first effect is the stan- 
dard energy decreasing due to the Universe's expansion; 
it is described by the term a 2 (to)/ a 2 (t) and tends to make 
the matter effectively nonrelativistic at t — > oo. A new 
effect caused by the nonstationarity of the dark energy is 
described in (|16[) by the exponential term. Let us empha- 
size, that this effect decreases the particle energy, if the 
term VV a ) [Il(t)— II(fo)] is positive, and increases it in the 
opposite case. Let us consider only one simple example, 
when V( a ) is positive, the initial value Tl(to) is vanishing: 
one can see, that the DM particle energy grows due to 
the exponential term, when the DE pressure is negative. 
In this sense, the Archimedean-type force with negative 
DE pressure counteracts the effective cooling of the dark 
matter in the course of expansion. In other words, even if 
a particle has small (but not equal identically to zero) ini- 
tial kinetic energy, it can become ultrarelativistic due to 
the Archimedean-type interaction with the dark energy. 
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3. Macroscopic moments of the distribution function 

Solving Eq. ([9]) taking into account (H"2"j) and (fT5|) we 

obtain 

/(»)(«, Pa) - / ( ° a) (g 2 )exp{3V (a) p(t)-n(t )]} , (19) 

where /(° a )(<7 2 ) is an initial distribution function of the 
DM particles of the sort (a). Using ([Ti?]) and (HT|) one can 
reduce the DM particle stress-energy tensor (the macro- 
scopic moment of the second order) 



ii DM) (t) = E/ 

(a) J 



dp 1 dp 2 dp 3 
V~9 Po 



f(a,)(t,P a ) PiPk (20) 



to the pair of basic integrals for the energy E and pressure 
P, respectively: 



In this model the DM energy and pressure have, respec- 
tively, the form 



pOO 

E ( X )=T,-1T q 2 dqJl+q 2 F (&) {x) e- A <*>V^, 

(a) X J ° 

(27) 

where the constant P( a ) is given by the formula 

_ iV( a ) 7Tl( a ) A( a ) 



(28) 



(29) 



The formulas for the Fermi-Dirac and Bose-Einstein func- 
tions can be obtained analogously, and we do not write 
them here. 



E^Yt^T 1 / q 2 dqf? a) (q 2 )Jl+q 2 F {a) (x) 
(a) X J ° 



(a) 



3x 3 



y/l+fFftix) 



(21) 



(22) 



Here for the sake of convenience we have introduced a 
new dimensionless variable x and an auxiliary function 
F( a )(x) defined as follows 

x=4ft, F (a) (a ; )=^exp{2V (a) [n(l)-n(a ; )]}. (23) 



a(t ) 



The initial moment t=to corresponds to the value x=l 
(F( a )(l)=l). The DE pressure can be also rewritten in 
these terms as 



n(x) = n(i) 



i 



2V, 



(a) 



log[x 2 F (a) (x)] 



(24) 



The most conventional model of matter deals with the 
relativistic Maxwell - Boltzmann functions [31j 



/(°a)(? 2 



iV (a) A (a) 



47rW (a)^2(A (a) ) 



-A (a )V 1 +9 2 



(25) 



where iV( a ) is the particle number density, A( a ) = k ™ ^ , 

and i^2(A( a )) is the modihed Bessel function [3l|, defined 
as 

/>oo 

2fn(A( a )) = / cosh vz ■ exp [— A( a ) cosh z] , (26) 
Jo 

k(B) is the Boltzmann constant. We assume that there 
is no thermodynamic equilibrium between different sorts 
of DM particles, thus, generally, the temperatures T( a ) 
do not coincide and are marked by the index of the sort. 



4- Massless particles mr a ) = 

The dark matter can consist of massless and massive 
particles. When we deal with massless particles, i.e., 
™(a)=0, e.g., for the index (a)=(0), we have to replace 
formally the quantity j?V TO (a) by K 1 with k l ki—Q. Respec- 
tively, the formula (|29|) has to be modified as follows: 



E, 



(0) 



47rt>[fc ( B)T , ( )]' 
h 3 



q 3 dq 
e« ± 1 



(30) 



Here h is the Planck constant, T(o) is an initial temper- 
ature of the massless particles, and v is a degeneracy 
factor. Clearly, this massless constituent of the DM is 
described by the ultrarclativistic equation of state 



E (0) (x) = 3P m (x)=E { 



Wo) 



(31) 



When the Archimedean- type force is absent, one obtains 
that F(q)(x)=x~ 2 and Em\{x) cx x~ 4 , as it should be. 

C. Balance equations 

1. Energy balance for interacting DE and DM 

We consider the model in which dark matter and dark 
energy form a coupled conserved system. The total 
stress-energy tensor of DM and DE is divergence-free 



rpih rpik I rpik 

-'(total) — J (DM) J (DE) 



V fc T (total) - ' 



(32) 



We suppose that both cosmic substrates, DM and DE, 
have the same macroscopic velocity four-vector, U l =5} jl 
thus, due to the FLRW space-time symmetry only one 
equation among (I32p is nontrivial: 

p + E + 3H(p + E + n + P) =0. (33) 

It is the direct differential consequence of Eqs. and 
©• 
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2. Balance equations for dark matter 

Balance equations for DM can be obtained by integra- 
tion of the kinetic equation (JT]) (see [3(| HJ for details) . 
The first-order macroscopic moment N k is defined as 



yp 



(34) 



and describes the total DM particle number four-vector. 
Using the kinetic equation (J7J the four-divergence of this 
vector can be easily calculated: 

V fc A fe = V fc ^ J dPp k f {a) = 



(a) 



d 

= -E m M J dP Qpk fi^k 
(a) 



= 0. 



(35) 



This means that the DM particle number is conserved. 
Taking into account the symmetry of the model one can 
write 



N k = S k N(t ) 



g(to) 
a(t) 



(36) 



The second-order macroscopic moment, the DM stress- 

pik 

-(DM)' 



energy tensor T?f$ M \ , satisfies the following equation 



(a) (a) 



v fc n^™^ a) v (a) f dPf {3 

(a) 



Ak P l P k 



P l Pl 



= S 3IT^V (a) F (a) . 

(a) 

Only one equation among Q37p is nontrivial: 
E + 3H{E + P) = -Q. 



(37) 



(38) 



The source term in the right-hand side of this equation 
has the form 



Q^3ri£v (a) P (a) 

(a) 



(39) 



The source Q in the DM energy balance equation van- 
ishes, when the DE pressure II is constant and the 
Archimedean-type force disappears. 

The scalar of the DM entropy production 

<T(DM)=V fc S fc = -fc (B )V fc ^ f dP p fe /(a) [l0gft 3 /(a)-l] 

(a) 



is also proportional to the four-gradient of the DE pres- 
sure 



C(DM)=fc(B) ' 
(a) 



'(a) 



=-3fc (B )V fc n^v, 



(a)%) 



(a) 



-3fc(B) n 



a(t ) 



1 3 



a(t) 



E V (a)^(a)(*0). (41) 

(a) 



When the coefficient 3? = J2(n) ^(a)-^(a) (*o) is positive 
and II > 0, the entropy production scalar <T(dm) is neg- 
ative, i.e., the DE provides organization processes in the 
DM system by the Archimedean-type force; when II < 
the Archimedean-type force produces a chaotization in 
the DM system. 



Balance equations for dark energy 



The combination of Eqs. (|3"3")l and (|3"5|) gives the fol- 
lowing balance equation for the DE state functions p and 
II: 



p + 3H{ P + n) 



(42) 



The functions P( a )(t) depend on II(i) according to Eqs. 
(f28l) and (j23|) for the massive particles, and according 
to Eq. (|3"Tj) for the massless ones. One can emphasize 
that the macroscopic balance equations (l38l) . (|39|) and 
(|42l) look like the balance equations in the well-known 
two-fluid models fl^ - Elj ]: the difference is that now the 
source term Q (j3"5|) is not modeled phenomenologically, 
but is directly calculated on the basis of kinetic approach. 



D. Dark energy dynamics 

1. Dark energy equation of state accounting for retardation 
of response 

To describe the dark energy fluid we use the linear 
equation of state 



p{t) = po + all + 



H(t) 



n. 



(43) 



It belongs to the class of the so-called inhomogeneous 
equations of state, which is intensely discussed in the 
literature (see, e.g., jH^). When cr=0 and £=0, Eq. 
(|43[) introduces the model in which the dark energy re- 
lates to the cosmological constant A [see, e.g., (d))]. Al- 
ternatively, po can be introduced by analogy with the 
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so-called bag constant appearing in the theory of quark- 
gluon plasma [38|]. When po=0 and £=0, Eq. (14*31 gives 
the well-known linear relation U=wp with w = — . Since 
the proportionality coefficient w depends on the choice of 
the epoch in the Universe's evolution, many authors con- 
sider it as a function of cosmological time, i.e., w = w(t), 
thus introducing the nonstationary equation of state (see, 
e.g., [32H37[ |). We follow another version of nonstation- 
ary equation of state, for which w and a remain constant, 
but the retardation of response is taken into account by 
inserting the term containing the first derivative of the 
pressure II. An equivalent scheme is widely used in the 
extended thermodynamics and rheology [see, e.g., [Hj]], 
in which the extended constitutive equation for the ther- 
modynamically coupled variables X and Y has the form 



rX + X = wY . 



(44) 



Here r is a relaxation time, a new coupling parameter of 
the model. In the cosmological context r is generally the 



function of time, r(t). We assume that r(i)= 



Tim- 



this relaxation time can be measured in natural cosmo- 
logical scale [4(j . Our ansatz here is that the dimension- 
less parameter £ is constant. 



2. Key equation of evolution of the DE pressure 

When the quantities pit) and II(t) depend on time 
through the scale factor a(t) only, i.e., p=p(a(t)), U — 
II(a(i)), the so-called x-representation is convenient, 
based on the following relations 



The quantity Jix) vanishes, when all the Archimedean 
parameters vanish, i.e., V( a ) = 0. Below we will address 
Eq. (j47l) as the key equation. It is a differential equation 
of the second-order linear in the derivatives but nonlinear 
in the unknown function II(x). There are two important 
particular cases, when the sourceterm Jix) can be writ- 
ten in an explicit form; let us consider them in more 
detail. 



3. Two explicit examples 
(i) Massless dark matter 

When one deals with massless particles [enumerated, e.g., 
by the index (0)], the source term Jix) takes explicit 
form 

Jp) (a:) = E m V(0) " {X) exp{V (0) [n(I)-n(x)]} (49) 



i.e., with Ei \ given by (|30[) . When DM particles are mas- 



(0) 

sive but effectively ultrarelativistic [q 2 F^ix) >> 1 and 
q 2 >> 1], the source term Jix) has the same form (|49"1) . 

(ii) Cold dark matter 

When one deals with the models using the concept of 
cold dark matter, one assumes that the corresponding 
particles are effectively nonrelativistic, i.e., q 2 F^ << 1 
and q 2 << 1. The corresponding key equation contains 
the source term 



= 3N {c) T {c) V (c) U'ix) e2V(c)[ri (i)-n(*)] _ 



(50) 



at ax 



t-tn 



'(*o) dx 



xH(x) 



(45) 



In these terms the balance equation (|4"2")l takes the follow- 
ing form (the prime denotes the derivative with respect 
to x) 

xp\x) +3i P + n) = 



(a)" 



[x 2 F {a) ix)]' q*dqe- x ^V^~ 2 



(a) 



2a; 4 



^l+q 2 F {a) ix) 



(46) 



Using (|43j) one can transform this equation into the equa- 
tion for the DE pressure only, yielding 

£x 2 n" ix)+xW (x) m+a) +3(l+(j)n+3/7o=7(x) , 

(47) 

where the source Jix)=Jix, II— 11(1), II') is defined as 
follows: 



(a) 



(48) 



Ultrarelativistic and nonrelativistic models can be stud- 
ied analytically and qualitatively. When the parameters 
A( a ) are arbitrary, one needs numerical analysis for the 
key equation. 



III. TWO EXAMPLES OF EXACT SOLUTIONS 

A. Constant solution to the key equation 

The behavior of the function n(x), the solution to the 
Eq. (|47p. essentially depends on the initial data 11(1), 
II'(l), on the values of the parameters A( a ), c, po, £ and 
V( a ). When a ^ — 1, the key equation admits a special 
constant exact solution 

n(x) = n(i) = -A = n(oo). (51) 

1+0" 

For this solution the energy density pix) is also constant 



P ix) = = -n(x) 



Po 

1 + CT 



-n(oo), (52) 



and the parameters ^, Xu\ and V( a ) are arbitrary. For 
this special constant solution the Archimedean-type force 
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vanishes (or more precisely becomes hidden) , and the cos- 
mological model turns into the FLRW-type model with 
dark matter and nonvanishing cosmological constant A. 
For this model we should write 



n(x) = - P (x) 



A 



and assume that 



Po 



A 



8nG 



(53) 



(54) 



taking into account the solution (fSlj) . This constant solu- 
tion also exists when the Archimedean-type force is ab- 
sent, i.e., when V( a )=0 and J{x)—Q. In addition, this 
solution is an asymptotic limit for a family of integral 
curves at x — > oo. 

A question arises: for what values of the parameters 
<7, po, £ and V( a ) is this special constant solution asymp- 
totically stable? In other words, when does the deviation 
from the asymptotic value tend to zero, Z(x — > oo) — » 0, 
where 



n(x)=- 



Po 
1+cr' 



-Z(x), Z(1)=0, Z'(1)=II'(1) ^0. (55) 



In order to answer this question let us analyze Eq. (|47[) 
linearized with respect to Z [see (|55[) ] at a; — > oo. In this 
case the leading-order term in the decomposition (|55j) 
satisfies the equation 



ix 2 z' 



(x) +xZ'(x) (4£ + cr) +3(1 + cj)Z = 0. (56) 



This equation is clearly the well-known Euler equation. 
First, it can be obtained as an exact solution for the 
case V( a ) = 0; second, this equation describes integral 
curves slightly deviating from the constant solution (fSTj) 
at x — > oo. The characteristic polynomial of the Euler 
equation (|56|) has two roots 



Sl,2 



1 

2i 



-(<t + 30± V(*-30 2 - 12£ 



(57) 



which can be real or complex depending on the values 
of the parameters £ and a. One can distinguish three 
subcases. 



1. Two different real roots [(a — 3£) 2 > 12£/ 
When the discriminant in (157[) is positive, one obtains 

(a; r - x~ r ) , (58) 



z(x) = m x -. 



2T 



where 



1 



7^^^, T^^(a-30 2 -m. (59) 



The solution (|58|) tends to zero asymptotically at x — > oo, 
when r < 7; it is possible, first, when a > 3£+2y / 3£ (for 
arbitrary £), second, when — 1 < a < 3£— 2 v / 3£ for £ > |. 



2. Double real roots [{a - 3£) 2 = 12£/ 

In this case the solution for Z(x) is 
Z(x) = n'(l) x" 7 loga;, 



(60) 



where 7=3+^/|, if er=3£+2V3£, and 7 =3-W|, if 
cr=3£— 2-\/3£. The first solution relates to the asymp- 



totically vanishing Z(x — > oo) for all £, the parameter cr 
being positive. The second solution tends to zero asymp- 
totically, when £ > i, and the parameter a satisfies the 
inequality a > — 1. Finally, when and <r=— 1, the 
solution is logarithmically unstable. 



3. Complex roots [(a - 3£) 2 < 12£7 
For the complex roots the solution is quasiperiodic 

Z{x) = 5111 x- 1 sin (ft log x) , (61) 

where 
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<t + 3£ 
2£ 



/3: 



1 

2? 



(62) 



This solution remains quasiperiodic and asymptotically 
small, when 7 > 0; it is possible in the following two 
cases: 

(i) when < £ < \ and -3£ < a < 3£ + 2^; 

(ii) when £ > i and 3£ - 2 v / 3£" < cr < 3£ -I- 2^. 

A special case, when the roots are purely imaginary, re- 
lates to the condition 7 = or equivalently a = — 3£. 
This special case is realized, when £ < -|. 

To conclude, we can state that the solution 
IL(x)—— is asymptotically stable, first, when < £ < 
I and a > — 3£, second, when £ > 3 and a > — 1. This 
statement is illustrated by the Fig.l. 

Since the zone of stability relates to the case a > — 1 , 
the asymptotic value of the DE pressure is negative and 
the dark energy is positive at po > 0. In the asymp- 
totic limit x — > 00 the energy density of the dark mat- 
ter asymptotically vanishes, E(oo) = 0, thus the Hubble 
function and the scale factor tend to the following values 



ff(oo) 



/ 8nGp 
3(1 + a) 



a(t) -> a(i )e H(oo) * 



(63) 



Such a behavior corresponds to the de Sitter law with 
positive acceleration parameter — q = — 1, thus, one 
can indicate this exact solution as a A - type solution of 
the model with Archimedean-type interaction between 
DM and DE. 

When g < — 1, the constant solution ([5~Tj) is unstable 
for arbitrary £. The case <r=— 1 is the special one, and 
we consider this case in detail in the next section. 
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FIG. 1. The stability domain of the special solution 
p=— II— presented on the plane of the parameters £ and 
a. The deviation from this constant solution decreases accord- 
ing to the power law (|58[) . first, in the zone situated higher 
than the upper branch of the parabola, a— 3^+2^/3^ for arbi- 
trary £, second, in the zone between the straight line a—— 1 
and the lower branch of the parabola, cr=3£— 2y / 3£ for £ > |. 
In the zone inside the parabola higher than the straight line 
o~= — 3£ the deviation from the constant solution behaves as a 
damped sinusoid (|6ip . On the branches of parabola the de- 
creasing is described by the product of a power-law function 
and logarithm (|60p . 



B. Special model with a=— 1: 
Exact solutions of the anti-Gaussian type 



For this specific model the key equation (147)) does not 
admit constant solutions if po =^ 0. In order to analyze 
the new situation we consider, first, massless dark mat- 
ter, then the model with cold dark matter, and discuss 
exact solutions of the key equation, which now become 
logarithmic, providing the scale factor to be of the anti- 
Gaussian type. 



1. Special model for a massless dark matter 

When <7=— 1 the key equation of the second order for 
massless dark matter 



£,x 2 U" {x)+xU l (x){A£,~l) 



n'(x)£ (0) ^l e v ( o)[n(i)-n(x)] +3iOo=0 (64) 



(0)^3 



can be reduced to the first-order equation 

& [4 + V ( o) W]^- = W[3£-l-E (0) V ( o) z] + 3p (65) 
based on the following definitions (see, e.g., fill]): 



dy 
dx 



y = IL(x) - 11(1) , x^- = W, z = x- 4 e- v ^ y . (66) 



Eq. (|65j) is the Abel equation of the second kind [41|. 
The special solution to this equation can be obtained if 
we put the coefficient [4+V(o) W] in front of the derivative 
equal to zero. Two constants 



W 



1 



V, 



(o) 



E, 



(0) 



3g-l 
V(o) 



;Po 



satisfy the key equation, if 



y(x)=-—— logx, 2=1, £(o) = 



(0) 



3£-l 3 



(67) 



(68) 



Of course, this special solution exists when V( ) ^ 0, i.e., 
when the Archimedean-type force is present. 



2. Logarithmic solution for the DE pressure 

When the parameters of the model are coupled by the 
last relation in f|68[) . we obtain an exact solution 



n(aO = n(i) 



v, 



(0) 



log X , 



p[x) = p(l) + -—logx 



V, 



(0) 



where the relation 



n'(i) = -^ = i[p(i)-p + n(i)] 

V (0 ) S, 



(69) 
(70) 

(71) 



is used. The sum IL(x)+p(x) remains constant for all 
values of x 



p(x)+U(x) = p(l)+U.{l)=po- 



4£ 



(0) 



The DM energy density E{x) is also constant, i.e. 

3 3^ i 

E(x) = const = E(q) = —~:Po ' 



V 



(0) 



(72) 



(73) 



due to the specific structure of the corresponding 
Archimedean-type force. The Hubble function H(x) can 
be found from Eq. ©, which now takes the form 



H\x) 



8ttG 



(0) 



V, 



logx 



(0) 



(74) 



Since the right-hand side of this equation has to be non- 
negative for every x > 1, we suppose that in addition to 
V(o) > the following inequality holds 



p(i) + e (0) >o ^n(i) < - Po 



v, 



(0) 



(75) 



According to (1431) and ([74]) the scale factor a(t) can be 
written in the form 



a(t) = a*exp<j ^-(t-t*) 2 



3V, 



(o) 



(76) 
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where the parameters with asterisks are defined as follows 



a* = a(to) exp 



V(o) 



p(l) + E (0) ] \ , 



t* = to — V( ) 



327TG- (0)J - (77) 

Such a scale factor describes an anti-Gaussian expansion 
without initial singularity. The acceleration parameter 
—q(t) for the anti-Gaussian expansion is positive and ex- 
ceeds one: 

3V (0 ) 



q(t) 



1 



> 1. 



aH 2 l6irG{t-t*) 2 

Moreover, it is a regular function of time (since t > to > 
t*), and tends to one, when t — > oo. Comparing the anti- 
Gaussian function (|76l) with the corresponding de Sitter 
type function 



(79) 



a(t) =a(t*)exp{\l-(t-t*) 



one can see that the former increases slowly at t ~ t* , 
but then grows much more quickly. 



3. Stability analysis of the anti-Gaussian solution 



The differential equation (|65[) can be evidently reduced 
to the autonomous dynamic system 

dz 

^=W[^-l-E {Q) V {Q) z] +3 Po . (81) 



tz[4 + V (0) W] 



(80) 



The exact constant solution (|67|) 

4 



IV* 



V(o)' 



2 = 1 . 



(82) 



valid at the condition (j68| , describes the stationary (sin- 
gular) point of this dynamic system. In order to deter- 
mine the type of this singular point let us consider the 
linearized system 

^ = £V (0) w , — = - Wj0o V (0 ) + 4£ (0) C , (83) 



where 



W^W„+uj, z^f z*+(. 



(84) 



Since we consider the case, when £ > 0, Er Q \ > 0, V( ) > 
0, the roots of the corresponding characteristic equation 
are real 



S l,2 



1± Wl 



256^ (0) 



9V ( o)P^ 



(85) 



The product of the roots Si ■ s 2 =— 4£-E(o) V(o) < is neg- 
ative, thus, the singular point is the saddle one, and the 
anti-Gaussian solution is unstable. 



4- Cold dark matter 



The key equation ([47]) with the nonrelativistic source 
(|ST)1) can also be reduced to the first-order equation 

Zz[5+2V ic) W}^-=W [3£-l-3iV (c) r (c) V (c) z] +3p 

(86) 

based on the following definitions: 

y=n(.x)-n(l) , x^-=W, z=x- 5 e- 2V ^ v . (87) 
ax 



ryg) The corresponding exact special solution is 



when 



G 



n(aO = n(i) 



2V, 



(C) 



■ log X , 



V 



(C) 



-(3£-l). 



(89) 



The corresponding value of the DE energy density is 



p(x) = p(l) + 



2V, 



(C) 



• log X , 



(90) 



so, that the sum of H(x) and p{x) again remains constant 
for all values of x 



P {x) + u{x) = P {\) + n(i) = po - 



5£ 



2V, 



(C) 



(91) 



The DM energy-density function E(x) is now a decreas- 
ing function of x (not a constant contrary to the massless 
case) 



E(x) 



m{c)N(c) 



(92) 



and the Hubble function H(x) can be now found from 
the equation 



H 2 {x) 



8irG 



2V, 



5 m iC )N {C ) 
log a; H 



(93) 



In the asymptotic regime x — > oo the logarithmic term in 
(1931) dominates, thus, we obtain again the anti-Gaussian 
solution 



aM = a**expj|^-n 2 



(94) 



where the parameters with double asterisks are defined 
as follows 

a = a(t ) exp i \ , 



t* 



to - 
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V(C) . /3p(l) 



2ttG 



(95) 



Again it will be an accelerated expansion of the Universe 
with —q > 1, and this solution is also unstable. 
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5. Model with DE domination 

In order to complete the analysis of the model with 
<j=— 1, let us consider the case, when the dark matter is 
absent (£V a )=0). In fact such a model can be considered 
as the approximate one, since at x — >■ oo the DM con- 
tribution to the total energy decreases as x~ 3 . The key 
equation for the DE pressure reduces now to the inho- 
mogeneous Euler equation 

£x 2 Il"(x) + xU'(x) (4£ - 1) + 3/?o = . (96) 

The characteristic equation for the corresponding homo- 
geneous Euler equation 



s(s — v) = , 



(97) 



with v = gives double roots s=0, when £=1/3, 

thus, let us consider two different cases. 
(i) Special case: a = — 1, £ ^ 

Exact solution to (|M|) contains in this case a sum of log- 
arithmic and power-law terms 



ncc)=n(i)- 



3j3 loga; (x u -l) 
l-3£ 1-3$ 



p(i)+n(i)- 



3£-l 



p(x)=p(l)-3^[n^)-n(l)]-3p logx. (98) 

Asymptotic behavior of the DE pressure and DE energy 
density at 3$ > 1 for arbitrary initial data are the follow- 
ing 



IL(x -> oo) — >• 



3|Qq 
3£-l 



log a;, 



p[x -> oo) -> _ - log a; . 



(99) 



In the asymptotic regime we deal again with the anti- 
Gaussian expansion 



a(t oo) oc exp i — _ t 



(100) 



when £ > i. We would like to mention that when the 
initial data 11(1), p(l) and n'(l) take special values 

n(l) = p(l)=0, n'(l) = — (101) 

the power-law terms disappear from the exact formulas 

11(00 = ^(1+3 log x), p(x) = ^Llagx, (102) 
and the anti-Gaussian solution 

a (t) = a (i )cxp||^L(t-t ) 2 | , (103) 



become not only asymptotic, but the exact one. When 
v > 0, then p(x) oc x v , H oc 2:2^ and a(i — > 00) — > 0. 
Thus, the Universe does not expand, when £ < |. 
(n) Special case: u=—\, £=3- 

The root s—0 of characteristic Eq.fjHT]) is now double, 
so, the exact solution contains the logarithmic terms in 
square, 

n(a;)=n(i)+3 [p(i)+n(i)- Po ] io g .T-^ log 2 x , (104) 



^ log 2 x, (105) 



p(x) = p(l)-3[p(l) + II(l)] logx 

the sum of these functions being linear in logarithm, 

IL(x) + p(x) = n(l) +p(l)- 3 Po log x . (106) 

The scale factor evolves now superexponentially. For 
instance, when the initial data satisfy the condition 
/9(l)+n(l)=0, the scale factor has the form 



a (to) 



= exp ■ 



' 2/7(1) 
9po 



sinh 



y/l2nGp (t-t Q ) }. (107) 



Near the starting point to the function a(t) behaves ac- 
cording to the de Sitter law 



a(t) * a(t )e H ^-^ , fffo) = ^ , (i 08 ) 
while asymptotically at t — ¥ 00 the scale factor grows as 

a(t) -> o(to) exp j e V^G P o (t-t ) | . (109) 

For the Universe expansion described by the law (|107[) 
the Hubble function is monotonic 



H(t) 



8nGp(l) 



cosh 



^/UnGpo (t - t ) (110) 



and increases infinitely. The acceleration parameter 
- <?(*) = 1 



9p Q sinh yUvGpo (*-*o)] 



2p(l) cosh 2 [V127rGpo (*-*o)] 



(111) 



starts with — q(to)=l, reaches 



"9(r. 



= 1- 



9po 



8p(l) 

asymptotically to — g(oo)=l 



at t ^o+ 1 °g (1+ ^ ) 

V 127rG Po 



the maximum 
and tends 



6. Remark on the coupling of the baryon matter to DE 

The baryon matter can be naturally included into the 
scheme of Archimedean- type coupling: for this purpose 
one can add to the sum the terms, which correspond to 
the standard particles. If one assumes that the baryon 
matter is not affected by the Archimedean-type coupling, 
one can put the corresponding coupling constants V( a ) 
equal to zero; nevertheless, the contribution of the baryon 
matter to the total stress-energy tensor now will be taken 
into account in the modified formula (|20[) . 
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IV. DISCUSSION 

The study of the cosmological model, into which the 
Archimedean-type interaction between dark energy and 
dark matter is introduced, shows that the roles of DM 
and DE in the energy balance of the Universe can be 
revised. According to the obtained formula (f!?Tj) , the 
contribution of the DM particles, E^(a(t)), into the 
total energy density depends on the state of DE pres- 
sure U(a(t)) through the functions F^(x) (see (f2"3"|) ). In 
the models without Archimedean-type force the energy 
of the DM particle decreases effectively because of the 
factor a~ 2 (t)\ in other words, all the particles, both non- 
relativistic and ultrarelativistic at the initial moment 
inevitably become (effectively) nonrelativistic in the pro- 
cess of Universe expansion. When the Archimedean-type 
force acts on the DM particles, the particle energy (fll))) 
becomes much more complicated function of cosmologi- 
cal time due to the exponential dependence on the DE 
pressure. This means, in particular, that, nonrelativistic 
particles can become (effectively) ultrarelativistic due to 
the Archimedean-type force action, thus the correspond- 
ing contribution of cold dark matter into the total energy 
can be reestimated taking into account the sign, the value 
of the DE pressure at this moment, as well as the rate 
of its variation with time. In contrast, the ultrarelativis- 
tic DM particles can become (effectively) nonrelativis- 
tic, when the corresponding exponential factor in (1161) is 
small. Now the contribution of cold dark matter into the 
total energy is estimated to be about 23%. The ques- 
tion arises: does this estimate include a total rest energy 
of the massive DM particles only, or the energy of the 
Archimedean-type interaction as well? We need such a 
clarification, for instance, in order to estimate the density 
numbers of the DM particles of different sorts; in partic- 
ular, the information about the number density of DM 
axions in the terrestrial laboratories is very important for 
plan ning experiments with axion electrodynamics (see, 
e.g., Ef). 

The cosmological model with Archimedean-type force 
describes a self-regulating Universe. This means that 
in the process of expansion of the Universe the total 
(conserved as a whole) energy can be redistributed be- 
tween dark energy and dark matter constituents accord- 
ing to the challenge of the corresponding epoch. The 
energy pendulum stimulated by the Archimedean-type 
force can work by the following scheme: let us imagine 
that at some moment the DE pressure II is negative and 
is varying rather quickly; then according to the formulas 
(1481) . (|49j) . ((50)) the DM particle reaction, provoked by 
the Archimedean-type force, will be strong. The corre- 
sponding intensive source appears in the right-hand side 
of the key equation (|4"T)) , thus decreasing the rate of IT 
evolution. From the theoretical point of view, it is not 
yet clear, first, for which set of guiding parameters such a 
specific regime does exist; second, when such a regime can 
be characterized as (quasi)oscillations; and third, how the 



number of epochs of the Universe expansion does depend 
on the set of guiding parameters of the model. We started 
to study these questions qualitatively and numerically in 
the second part of our work and presented examples of 
quasi-periodic, multi-inflationary and multistage evolu- 
tionary schemes in the framework of the model based on 
the Archimedean-type interaction between dark energy 
and dark matter. 

In the first part of the work we focused on exact so- 
lutions of this new model. The first exact solution is 
the constant one, U(x) =— p(x)=— and relates to the 
case c=^—l. Since the DE pressure for this exact so- 
lution is constant, the Archimedean- type force becomes 
hidden, and we obtain the standard cosmology with A 
term. This solution is asymptotically stable, when the 
guiding parameters of the model, £ and <r, satisfy the in- 
equalities < £ < | and a > — 3£, or £ > | and a > — 1. 
When a < — 1, the solution is asymptotically unstable 
for arbitrary parameter £. 

Exact solutions of the second class, with a = — 1, 
are much more sophisticated, since the DE pressure 
is described by the logarithmic function of the ratio 
a(t)/a(to). The corresponding scale factor a(t) is given 
by the anti-Gaussian function (|76p . it has no singular 
points in the early Universe, and describes late-time 
accelerated expansion with the acceleration parameter 
— q > 1; this acceleration parameter is bigger than for 
the de Sitter model. Exact solutions of the anti-Gaussian 
type happen to be typical for different physical situations: 
for massless and massive nonrelativistic DM, for the case 
with DE domination, etc. This solution is unstable. 

There are two specific values of the guiding param- 
eters of the model: <r=— 1 and The first one, 
o=— 1, can be associated with the so-called phantom di- 
vider w(t)=— 1=^. The second value, £=^, can be de- 
noted as a resonance value of the relaxation time param- 
eter. Indeed, according to the formula (|43|) the func- 
tion -gl^j plays a role of the relaxation time for the DE 
pressure TI(t), say, tu- When £=^, one obtains that 
Tn= W(t) ~ep) ' wnere Q(t)=3H(t)=VkU k is the expan- 
sion parameter. Thus, the characteristic time of expan- 
sion coincides with the relaxation time parameter 
for the DE pressure, introducing some specific resonance 
condition. When a—— 1 and £=i simultaneously, there 
exists superexponential solution of the model, described 
by the scale factor (|107D . 

A physical origin of the Archimedean-type force is not 
yet clear; nevertheless, this force seems to be very inter- 
esting from the viewpoints of a new model of interaction 
and a new scheme of redistribution of the cosmic energy 
between the interacting DE and DM constituents. The 
presented model is self-consistent, simple from the point 
of view of analysis and very promising. We keep in mind 
the story of the Chaplygin gas model 43] : starting from 
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a classical analogy a new evolutionary model has been 
elaborated and applied to cosmology, although physical 
meaning of the Chaplygin scheme of interaction is under 
discussion till now. 
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